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Abstract 

We present a calculation of the two-loop anomalous dimension for the transversity dis- 
tribution hi{x,Q'^), 7^*-^'', in the MS scheme of the dimensional regularization. Due to the 
chiral-odd nature, hi does not mix with the gluon distributions, and thus our result is the 
same for the flavor-singlet and nonsinglet distributions. At small n (moment of hi), 7^*-^'' 
is significantly larger than 7i{*^^^ (the anomalous dimension for the nonsinglet /i), but ap- 
proaches 7^*-^^ very quickly at large n, keeping the relation 7^*-^'' > 'jI^^^ ■ This feature is in 
parallel to the relation between the one-loop anomalous dimension for /i and hi. 



PACS numbers: 13.85.Qk, 12.38.Bx, 13.88. +e, 14.20.Dh 



1 Introduction 



It is well known that the nucleon (spin-1/2 hadron) has three independent twist-2 quark 
distributions fi{x,Q'^), gi{x,Q'^) and hi{x,Q'^). These distributions have a simple parton 
model interpretation and appear as a leading contribution to various hard (semi-) inclusive 
processes. The first two distributions, /i and gi, are chiral-even and can be measured through 
the totally inclusive deep inelastic lepton-nucleon scattering (DIS). In fact, accumulated DIS 
data on fi and gi has provided us with a rich information about the parton distributions in 
the nucleon. The third one, hi, is chiral-odd and cannot be measured through the totally 
inclusive DIS. However, it appears as a leading contribution to the transversely polarized 
nucleon-nucleon Drell-Yan process |^, 0, Q , semi- inclusive productions of baryons and 
pion in the DIS [|, |^ , etc. This spin-dependent quark distribution, named as transversity 
distribution in [Q, measures the probability in a transversely polarized nucleon to find a 
quark polarized parallel to the nucleon spin minus the probability to find it oppositely 
polarized. It thus supplies an information on the nucleon's spin structure not provided by 
the gi distribution, and is expected to be measured by RHIC in the near future. 

The Q^-evolution of these parton distribution functions is calculable in perturbative QCD. 
It is governed by the anomalous dimensions for the corresponding twist-2 operators and the 
/9-function. Leading order (LO) Q^-evolution of /i and gi has been known for a long time 0] 
and the LO Q^-evolution for hi was also studied in 0, For the next-to-leading order 
(NLO) evolution, we need two-loop anomalous dimensions. So far two-loop anomalous 
dimension was calculated for fi ||, ||, [ll], [12| and for gi [|l^, [l^. They have been 



used for the phenomenological parametrization of fi ||TB[ and gi [|TB| to analyze accumulated 
experimental data. 

The purpose of this paper is to calculate the two- loop anomalous dimension for hi, 
following the method of 0, |[. We calculate the two-loop correction to the two-point Green 
function which imbeds the twist-2 operators for hi{x,Q'^). The calculation is done in the 
MS scheme in the dimensional regularization. 

The construction of this paper is the following: In section 2, we briefly recall the twist-2 
operators for the distribution hi and summarize the formal aspect of the next-to-leading 
order (NLO) effect in the Q^-evolution of hi. In section 3, we describe the actual procedure 
of our calculation. Although the method is in parallel with that of |^, presence of the 
free Lorentz index /i in the twist-2 operators for hi brings additional complication. We shall 
describe this peculiar feature for hi in great detail, but omit the details which are common 
to both hi and fi. In section 4, we present the results and discuss characteristic features 
of the NLO Q^-evolution of hi{x,Q'^). Readers not interested in the technical details can 
directly proceed to this section. Section 5 is devoted to a brief summary. In Appendix A, 
we present the necessary integral formulas to express the anomalous dimension contribution 
from each diagram. In Appendix B, we give the contribution to the anomalous dimension 
from each two-loop diagram. 
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2 Twist-two operator for hi{x^Q'^) and its renormaliza- 
tion 

The transversity distribution hi{x, Q"^) of the nucleon is defined by the following relation 0: 

+hL{x,Q'^)M{p^n^ - Pi,nf,){S ■ n) + h3{x,Q'^)M{S_L^ni, - S±un^) , (2.1) 

where \PS) is the nucleon (mass M) state specified by the four momentum P and the spin 
vector S, and two light-like vectors p and n are introduced by the relation, P^ = p^ + 
and pn = 1. is decomposed as = {S ■ n)p^ + {S ■ p)n^ + S'^. In (|2.1|), the gauge 
link operator to ensure the gauge invariance is suppressed, and /is are twist-3 and -4 
distributions which are not in our interest here. Taylor expanding the both sides of (|2.1| ) 
with respect to A, one can relates the ra-th moments of hi ^ to the following local operator: 

Of,u|l^...t^,, = Sni'il->(r^,uiD^^ . . . iD^^i), (2.2) 

where Sn denotes the symmetrization among u, fii, ■ ■ ■ , fin, and the covariant derivative = 
— igT°'A1^ restores the gauge invariance. In particular, the n-th moment of hi is related 
to a tower of the twist-2 operators: 

J^^ dxx''hi{x,Q^) = a„{Q^) (2.3) 
(P^IVm...,„(Q)A^A^^---A'^"|P^) 

= {S,P'^+' - P^^P" + -l^M^A^SP"-!) , (2.4) 

where 6^^^^...^^ is defined as the traceless part of 6fj,u^^--fj.„ defined by the condition, 

Q^^""^ fj.vij.\...^in — 9 '^'^fiu^ii...^i„ — 9^ ^ iJ,u^ii...^i„ — 0, (2-5) 

and we introduced a null vector A'^ (A^ = 0) which kills off the trace terms in 6^^^^...^^ as 
usual. 

For simplicity we symbolically use the notation 0„ for in what follows in this 

section. The bare-(0;f) and the renormalized- (On(/i)) composite operators are related by 
the renormalization constant Zn{fi) for On as 

0„(/i) = Z-\fi)0^. (2.6) 

The scale dependence of Oni^jj) is obtained by solving the renormalization group equation 

/^^^ + 7n(^(/x))O„(/x) = 0, (2.7) 

where 'jn{9(yfJ')) is the anomalous dimension for the operator On^jj) defined as 

d 

7n.(^(/^)) = /^^lnZ„(/i). (2.8) 



This equation is solved to give 



On{Q^) = 0„(/i2)exp 



(2.9) 



Up to next-to-leading order, the anomalous dimension 7„(/i) and the beta function [3{g) can 
be expanded as 
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(2.10) 
(2.11) 



where the coefficients of the /^-functions are well known; j3o = 11 — {2/3)Nf, f3i = 102 — 
(38/3)A^/, with the number of quark flavor Nf. Inserting these equations into (p.9|) , one 
obtains the next-to-leading order evolution equation for the n-th moment of hi as 



1/ 7i^_2^ 

/30 



with 



9'm 



,/5o + /?i^?2(/i')/16vr2 



/?ilnln(gVA2 



(2.12) 



167r2 /?oln(QVA2) 



+ 



(2.13) 



where we introduced the shorthand notation Ain[hi{fi)] = j\ dx x^hi{x, This NLO 
effect in hi has to be combined with the NLO effect in the hard cross section in the parton 
level to give a prediction for a relevant physical cross section in the NLO level. 
In the MS scheme of the dimensional regularization, one obtains Zn{n) as 



= 1 + E 



(2.14) 



fe=l 



with e = {A — d)/2. Then the anomalous dimension is obtained as 



7„(^?) = -2g' 



(2.15) 



Therefore 7^^-' in ( p.lO| ) is —4 times the coefficient of the g^ / {Wn'^y term in X"(/i). (See 
for the detail.) 



3 Calculation of the two-loop anomalous dimension 

To obtain the renormalization constants for 0i^'^t^^" i^" in (|2.4| ), we calculate the two-loop 
correction to the two-point function for 6 with off-shell quark lines having the momentum p 
{p^ < 0). The calculation of the two-loop Feynman diagram is rather involved, and we follow 
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the procedure of [0, To discuss renormalization of 9, the presence of ^75 is irrelevant, since 
ij^cT^'^ = —l/2e^^^'^Pai^p. We shall henceforth omit for notational simplicity and consider 
the renormalization of the operator Oi^^^^'"^" defined by 



(3.1) 



where iST„ denotes the symmetrization among i/, /xi, ■ ■ ■ /i„ and the traceless condition in 
02.51). As usual, we contract Q^'^^^'"^'^ with a null vector A'^ to kill off the trace terms. 



and consider the renormalization of • A = O^'^'^^ "'^" Aj^A^^ ■ ■ ■ A^^ in the following. The 
calculation is done in the Feynman gauge for the gluon propagator, but the results should 
be independent of the gauge. 

We first need the Feynman rule for the basic vertices for the operator ■ A with 0, 1, 
2 gluon lines shown in Fig.l: 



(9(i)m . 



(3.2) 
(3.3) 



l<j<l<n 



l<j<l<n 



(3.4) 



where we introduced the notation p = p ■ A for a four vector p. If one replaces a^'^A^ by 
^ in (|3.2|) - (|3.4|) , one obtains the Feynman rule for the same basic vertices for /i, which are 
used in [0. Relevant two-loop diagrams are shown in Figs. 2, 3 and 4. Diagrams shown in 
Fig. 3 are only for the flavor- singlet distribution, but they are identically zero for hi owing 
to its chiral-odd nature. Those in Fig. 4 are the self energy corrections and are given in |^ . 

Among the contribution from each diagram in Fig. 2, we can show that the figs, (d), 
(g), (h), (1), (m), (q) and (r) give the same anomalous dimension to hi as the fi distribution 
calculated in H7|. To see this, we note the relation 



a^^'A^ = i^^^ - %A^ 



lA^ 



(3.5) 



The term proportional to A'^ is the trace term, and thus we can drop them and use either 
^7^^ or —i^'-ff^ instead of a^^Ay in the actual calculation. We shall henceforth use the 
three forms for the vertex factor, a^^'Aj,, 7^^^, —^7^, interchangeably (we also drop factor i 
for simplicity) without mentioning explicitly. Therefore the difference in the Feynman rule 
between hi and fi is merely the presence of i'^^ next to ^. The loop corrections in the above 
diagrams are located in the right or left half of the diagram and hence the presence of 7^ 
does not bring any change in the calculation of the renormalization constants compared with 
the /i case. We can thus take the results from for these diagrams. 

It is also easy to confirm that the diagrams (a), (b) and (c) become finite due to the 
form of the vertex, and hence no contribution to the anomalous dimension unlike the case 
for fi. Diagram (p) has no contribution as in the case of fi. Calculation of (e) and (f) is 
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also standard. One gets contribution from these two diagrams which are independent of n 
unhke /i case. 

Rather heavy calculation is required for the diagrams (i), (j), (k), (n) and (o). Compared 
with the renormalization of /i, additional complication is caused by the presence of the free 
Lorentz index /i in ■ A. To explain our strategy of the calculation, we take the diagram 
Fig. 2 (i) as an example. We assign the momentum as shown in Fig. 5. Application of the 
Feynman rule results in the following contribution from the diagram. 

D{k^, k^) = k{kl{k^ - p)\k2 - v)\k^ - k2)\ (3.7) 

iV(fci, k2) = jxh^^^'^uhipih - hh^i^ - (3.8) 

where 

(3.9) 

Our task is to extract from (|3.6|) the 1/e pole terms proportional to a'^'^A^p"'. To this end 
we first contract the index p in the numerator ( p.8| ), noting the properties of the 7 matrix 
in d = 4 — 2e dimension. This gives 

-Ae{k^ - k^) ■ {p - k2)ixh^^^A4^^\ (3.10) 

Noting 7A^icr'"'Ai,^i7'^ = -2elj^ia^"' A^J/^i, the last term of (|3.10|) becomes 0{e'^) and cannot 
produce 1/e-poles so that we can drop it. Using the same relation in the first term, one gets 
up to 0{e) 

N{k,, k2) = -4(^ - Iit^Wi^^^'^Mh - h) - 4£^icr^"A,^i(2i - ^2)(^i - h) 

+Ae{lir - mi^^'MS - h)- (3.11) 

In the first and the third terms in ( p.ll| ), we can discard the terms which have at the 
leftmost or the rightmost places, since those terms are either ultraviolet finite from power 
counting or only cause trace terms which are proportional to A'^. In fact, for example, one 
can write 

+5yp"+i(p2)-(i+2^) + CA^^^p'^-'ip^)-^' + C"A''p"(/)-2-, (3.12) 
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where A, B and B' terms are finite, and C and C terms are the trace terms. Dropping these 
terms, one can rewrite as follows: 



+Ae{2ki ■k2-2p-ki + 2p-k2- A;^)^i7''^^i - AeklljtiYM2- (3.13) 

The integration of the first and the second terms in ( |3.13D can be handled as follows. One 
can write 

J J D{ki,k2) 

where A can contain a pole at e = 0. The A term can be projected out by multiplying 
to both sides of ( p.l4| ) and taking trace of both sides. This is equivalent to replace in the 
numerator of (|3.14| ) as 

^2h^l'' ^{pki ■ k2 + hp ■ k2 - hp ■ ki)^-i^'. (3.15) 



Similarly to ( |3.15| ), we can replace in the second term of ( p.l3| ) as 



7"A^i|^2 ^ \{pki ■ k2 + hp ■ k2 - hp ■ ki)YA (3.16) 
P 



After the replacement (|3.15| ) in (|3.14|) , the integration can be easily done by the same method 



as For example, ki ■ ^2 in the first term of ( |3.15| ) is replaced by l/2{kl + k^ — {ki — k 



"2 J 



^2^ 



Then each of these three terms combined with the kf factor in ( |3.14| ) cancel some of the 
factors in the denominator D{ki, ^2) in (|3.7| ), leaving the 0(A;f 2) factor in the denominator. 
After this cancellation, one needs only one Feynman parameter for each integration (/ Dki 
and / -0^2) and thus one can easily carry out the integration. 

As is shown in this example, our strategy of the calculation is to devise a suitable replace- 
ment for each term in the numerator to project out the pole terms proportional to 7^^^ or 
^7^^, ignoring the trace term (~ A'^). After this projection, it turns out that we can express 
the final formula for the pole contributions of the two-loop diagrams in terms of the integral 
formula Ii to lu listed in the appendix of [Q and two new integrals /15 and /ig shown in the 
Appendix. 

As was stated in ( p.l5| ), simple pole terms (~ 1/e) from each two-loop Feynman diagram 
give rise to the anomalous dimension. There are two origins for the simple poles: (i) 0{e^) 
contribution in the numerator combined with the 1/e singularity from the integral, (ii) 0{e) 
contribution in the numerator combined with the 1/e^ singularity from the integral. For 
the type-(i) singularity, one can subtract the counter terms after transforming each two- 
loop integral into the form of Ji ~ Jig. Therefore, in the integral formulas given in ||^ and 
Appendix A, the counter terms for the simple pole terms are already subtracted. To subtract 
the counter terms for the type-(ii) 1/e-pole, we have to come back to the original diagram |^. 
In Fig. 6, we schematically showed the counter term contribution for the diagram (i) in Fig. 
2. In Fig. 6(a), the counter-term corresponds to the renormalization of the operator in 
one-loop (Fig. 6(b)), but the numerator of the diagram in Fig. 6(a) has no e. In Fig. 6(c) 
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the counter-term comes from the one-loop contribution to the vertex correction (Fig. 6(d)). 
This time the numerator of the diagram in Fig. 6(c) becomes 0{e^), hence no counter-term 
contribution. To summarize: we only have to subtract the counter-terms for the simple poles 
in Ji ~ Ji6 for the diagram (i). This is also true for other diagrams in Fig. 2. 

The integration of the third term in ( 3.13|) is rather demanding. We repeatedly apply 



the anti-commutation relation to move the leftmost ^2 untill it merges the righmost to 
form ^2^2 = ^2- C)ne then obtains 

The first term of ( p. 171) can be handled in the same way as (|3.16| ). 

The fifth term of ( p.l7| ) is the easiest to deal with. One only has to make the replacement 



7^^2^^7^A (3.18) 



by the same argument as in ( 3.14| ) 



The third term of ( |3.17| ) can be treated as follows. To meet with more general situation 
we replace /c2 by /cf (fcf can be put equal to fcf or k!^ in the following formula.) From power 
counting, one can write 

r r k'^k'^k^hUn 

/ Dk^ / Dfc2 i./, T = + (finite as e ^ 0), (3.19) 

J J D{ki,k2) 

where A and B can contain poles at e = 0. To project out the A-term, we multiply ^^7^ from 
the left and take the trace. To 0{e), this is equivalent to the replacement, 

k^M2 ^ ^(1 + £){kiP ■ k2 + hp ■ ki - pk, ■ k2)^-i^, (3.20) 



in the numerator of (|3.19| ). Similarly to (|3.20|) , we often meet the following replacement in 
the calculation: 

fc^^2A ^ ^(1 + e){kiP ■ k2 + hp ■ ki - pfci ■ A;2)7^A- (3-21) 
2p 



We can deal with the second and fourth term of ( p.l7| ) as follows. Without loss of 
generality one can write {q'^ = k^ or k^) 

/ Dki / Dk2 , ; = A7^A + BA^ + CA^j^A + (finite as e ^ 0), (3.22) 

J J D[ki,k2) 

where A, B and C can contain poles at e = 0, and we can ignore finite contribution. To 
project out the A-term, we combine two operations, (i) Multiply p^'/i^ from the right and 
take trace, (ii) Multiply 'y^^ from the right and take trace. Taking the difference between 
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these two results, it turns out that the projection of the A-term in ( p.22| ) is equivalent to 
the following replacement in the integrand: 

q%h ^ 7^(1 + ^) {PiP ■k2q-ki-p- hq ■ fca) - P ■ q{kiP ' h - kp ■ kMj^^. (3.23) 
Zp 

The sixth term in ( ^Tft j can be rewritten as ^iYM2 = 2A;f^^2 - 2kiYh + l^Mih and 
we can apply ( p.20D , (|3.18| ), (|3.16|) for each term, and likewise for the seventh term in (|3.17|) . 
As an application of the above replacement rules, we prove a usefull lemma here. 

Lemma: A term of the form eJ/^a^^ A^J/^' , where k and k' are either ki or /c2, does not give 
rise 1/e pole. 

Proof: We write 5^7^ A^' = 2ek>'^J^' - 2ekYf + e^i^W, and apply (pD, (|3A8|) , (|3A6| ) 
for each term. Then it vanishes, ignoring 0{e^) term. This completes the proof. 

Applying this lemma, we can drop fourth and the fifth terms in ( p.l3| ), which completes 
transformation of all terms in ( |3.13| ) into the requested forms as mentioned before. By this 
lemma we can drop many of the terms in the calculation of the diagrams (i), (j) and (k). 

The obtained rules for the replacement ( CT) , ( CT ), (|3l8| ), (lOOl) , (p%D , j ^M) give 



all the necessary formula to carry out the integration of the two-loop diagrams. Using the 
method described in this section, we obtained the anomalous dimension from all the diagrams 
in Fig. 2. They are given in Appendix B. 



4 Results and discussions 

Collecting all the contribution from each diagram in Figs. 2 and 4 (See Appendix B and 
appendices of 0), we obtained 7^^) for hi as. 



In 



where 



4C| 



^2(n+l)-25i(n + l)-- 



+CfCg 



^8 ^77 90 

-16^1 (n + l)S2{n + 1) - -S2{n + 1) + —Si{n + 1) - y 



8 I Cp — —CpCc 



ASi{n + l)\S: 



, fn + 1 



S2(n + 1) 



8S{n + l) 



, fn + 1 



-S' 



y 



-S2(n + 1) + \ + - 

2 ' n + 1 n + 2 4 



352(n + l)-5^i(n + l) + - 



Sk{n) 



Sh. 



n ^ 

n 



-11 



-Sh 



-Sh 



n — 1 



(4.1) 

(4.2) 
(4.3) 
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S{n) = E^^^i(j)- (4.4) 

j=i J 

(Note that 2/35'i(n) in the seccond hne of (A. 8) in should read 2/3S'3(n).) This final 
formula for hi{x, Q^) is very complicated as in the case of fi{x, Q^) and is subject to numerical 
analysis to study its Q^-evolution. In table 1, we present the actual numbers of 7^*^^^ (7^^^ 
for hi) for n = 0, ...,20. This table can be compared with Table 2 of 0, which presents 
7n''^'* i'ln^ /i)- (Note that n in [|^ corresponds to our n — 1.) We also plotted 'j^''^^^^ 
for Nf = 3, 5 in Fig. 7 to compare their behavior. From Fig. 7, one sees clearly that at 
small n 7^^"^-' is significantly larger than 'jI^^^ but approaches very quickly to 7^^^-', keeping 
the condition 7^^^-' > 'J^^^- This feature is the same as the one-loop anomalous dimensions 
/y^,/(o) -y^^iiich reads 

= 2cJl + 4E- , (4.6) 



i=2 



J 



and hence 7^*-°^ > 7^*-°'' for all n. Actually this tendency is even stronger for the two-loop 
case. We note that for n = the anomalous dimension for /i is zero in all orders because 
of the Ward identity for the vector current. On the other hand, hi projects onto the tensor 
operator ipa^^ip for n = 0, for which there is no conservation law. Hence 7o*'°'^'' 7^ 0. This 
can be taken as a good reason for 7^^°'^^ > 7i{^°'^'* especially at small n. From the difference 
between 7;(^^'' and 7^^^'' at small n, we expect that the NLO effect leads to larger difference 
in the Q^-evolution in the small-x region. 

At large n, 7^'-^'-^'* receives dominant contributions from diagrams (d), (g), (h), (m), (q), 
(r) in Fig. 2 which give the same anomalous dimension both for hi and /i, and sub-dominant 
contribution from (n), (1), (f), (e) and the self energy corrections shown in Fig. 4. All other 
diagrams in Fig. 2 gives minor corrections, although they are very important at small n. As 
is discussed in 7^*-^^ behaves as ~ ln(n) at large n with the cancellation among ~ In^n 
and ~ In^n terms which arise from some of the diagrams. Equation ( [4.1|) and quick merging 
between 7^*-^'' and 7^*-^'' at large n ( n > 10 ) with the relation 7^*-^'' > 7^*-^'' for all n clearly 
shows 7^*-^'' ~ ln(n), which also supports the correctness of our calculation. 

The relevant quantities for the Q^-evolution of the moments are 7^°Y2/3o and 7^^V2/3i 
as is seen from ( |2.12| ). In Fig. 8(a), we plotted these quantities for /i and hi with Nf = 3. 



Since the NLO effect is determined by their difference (see ( |2.12|) ), we plotted ^^'^^^y2(3i 
'jf^Mo) ^2Po in Fig. 8(b) for Nf = 3 and 5 cases. From Fig. 8(b), one expects that at small n 
the NLO effect in the evolution is quite different between Mn[hi{Q'^)] and -Mn[/i(Q^)]. 
Since 'jl^^'^^ — as n — 0, 7^*^^V2/?i ~ 7n^°V2/3o abruptly drops to zero as n — >• 0. But 
this is not the case for hi. The difference between the black circles and squares in Fig. 8(a) 
shows a characteristic behavior as shown in Fig. 8(b). 

As an example of the Q^-evolution, we show in Fig. 9 (a) and (b) the Q^-evolution of the 
tensor charge and the first moments, respectively, with the parameters Nf = 3 and A = 0.232 



GeV in ( 2.13 ). (Here we are interested in the NLO effect in the anomalous dimension and 
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the /5-function, and thus we adopted the same value for the A-parameter in the LO and NLO 
evolution.) They are normalized at = 1 GeV^. At n = 0, only diagrams (e), (f), (i), (j), 
(k) survive. They give the anomalous dimension for the tensor charge as 

= Nf. (4.7) 

9 27 ^ ^ 

From Fig. 9(a), we can compare the LO and NLO Q^-evolution of the tensor charge of 

the nucleon. One sees that the NLO effect is sizable as is expected from Fig. 8 (b). In 

Fig. 9(b), we plotted the LO and NLO Q^.evolution both for Mi[fi{Q'^)] and Mi[hx{Q'^)]. 

Although the NLO effect in the anomalous dimension (the second factor in the right hand 

side of ( |2.12| )) makes this ratio smaller, the NLO effect in the coupling constant (the first 

factor in the right hand side of (|2.12 )) completely cancels this effect. For /i, the latter effect 

is actually bigger than the former effect. As was stated in section 2, the NLO effect in the 

distribution function has to be combined with the NLO effect in the hard cross section in 

the parton level to give a prediction for a physical quantity. We will pursue this issue in a 

future publication. 



Table 1 The coefficients of the two-loop anomalous dimension 7,^^^ for hi with the number 
of quark flavors Nf = 3, 4, 5. 
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165 


79 


11 
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16 


248 


96 
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79 
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77 


224.55 
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32 
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258 


36 


228.53 


198 


71 
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232.32 


201 


92 
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266 


91 


235.96 


205 


01 
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5 Summary 



In this paper, we have carried out the calculation of the two-loop anomalous dimension for 
the transversity distribution hi{x,Q'^). The calculation was done in the MS scheme of the 
dimensional regularization. This completes the calculation of the anomalous dimensions for 
all the twist-2 distributions of the nucleon in the NLO level. We found 7^^^-' is significantly 
larger than 7^^^-' at small n, but approaches very quickly to 'y^^^^ at large ra, keeping the 
condition 7^^^-' > 7i{''^^- This means that the NLO effect for hi is quite different from that 
for /i in the small x-region. As an example, we have compared the LO and the NLO Q^- 
evolution of the tensor charge and the first moment of /i and hi. We hope the peculiar 
feature of the Q^-evolution of hi studied in this paper will be measured in the ongoing and 
planned experiments in the near future. 



Note added: 

After submission of this work for publication, there appeared a preprint [1^ which calcu- 
lated the NLO splitting function for hi by a different method. Our whole result agrees with 

which calculated 



p!7| . After submission, we also learned that there appeared a preprint 
the two-loop anomalous dimension for hi in the same formalism as ours, but presented the 
result in a very complicated non-final form. We further noticed that (1) their result for the 
diagrams Fig. 2 (i), (j), (k) disagrees with ours, (2) their result for Fig. 2 (n) and (o) agrees 
with our intermediate result but are not presented in a simple form. We believe something 
must be wrong in their calculation. 



11 



APPENDIX 



A. Integral formula 

To write down the two- loop anomalous dimension for hi, it turns out that we need two more 
integral formula, /15 and Jig, in addition to Ji ~ /14 defined in the appendix of [§. As in 
the case of Ji ~ /14, we define them by subtracting the counter term contribution for the 
simple-pole term but not subtracting the counter-term for the double-pole term as explained 
in the text: 



/i5(M,0) 

/i5(0,0) 

/16 (M,iV) 
Iw{M,N) 



kfki-k2 



{27r)D (27r)^ kf{ki - pfik^ - pfih - ^2)^ 

M 3 



p 



(167r2)2 4(M+ 1) 



1 1 
2^2 + 2i 1 3M ' M + 1 



1 1 

+ 



- 2 



-1 



4(167r 



2^2 



(M > 1) (A.l) 
(A.2) 



d^ki d^k2 



kfk^ki-p 



p 



{27i)D (27r)^ ki{k2 - pf{ki - k2Y{ki -k2+ pY 



E 



2(167r2)2 (M + A^)! (A^-A;)!(M + A; + 1)! 
1 



1 

2^2 



— { Si{M + N -k-l) + Si{M + k + l)- 2Si{M + A^ + 1; 



(M > 1, N>0) 



-p 



N 



N 



(47r)M(A^+ 1) 



1 1 



N-l 



k=0 



-2Si{k + l) + Si{k)-S_ 



N-k-l 



k + V 

k + 2 



{N>1) 



/i6(0,0) = 0, 



(A.3) 



(A.4) 



(A.5) 



where we discarded the finite contributions. We also remind the readers that our e is e/2 in 



Although the expression for Ji ~ /14 are given in P], some of them contain misprints. 
Furthermore, when we calculate the anomalous dimension for the tensor charge, we need to 
prepare some of Ii ~ /14 separately as in the case of /is.ie- In the following, we present such 
integral formulas. 



13(0, AT) 
h{M,0) 



p 



2(167r2)2 1 
d^ki d^k2 



\--SiiN + l] 



{N > 0) 



(A.6) 



k^k2 ■ p 



(27r)^ {2n)D kf{ki - pY{k2 - pf^ki - /ta)^ 



12 



1 



IsiM,N) 
h{M,N) 



4(167r2)2M + 1 



1 



1 f 

+ £ 1 2(M + 1) 



+ 



2M 



- 1 



(M > 1)(A.7) 



(27r)^ {2n)D kikl{k2 - V?{h - k2) 



P 



M+N 



4(167r2)2 {M+1){M + N) 
1 12 

+ 



1 1 
^ + - 



h{0,N) 
hi{M,N) 



M + N + 1 M 
-p^ 1 1 



d^ki d^k2 



M + 1. 

{N > 0) 



iMtN 
'^1 '^2 



^i(M + iV) - Si{M) 
(M > 1, AT > 0) 



(A.8) 
(A.9) 



{27i)D {2tx)D kl{k2 - py{k, - k2y{ki -k2+ pf 



-p 



M+N 



(167r2)2 e'^M + N + 2 



1 {-l)'^M\N\ 
M+1 ^ (M + Ar + 1)! 



(A.IO) 



The summation in Ji and /13 should start from J = in (A. 2) and (A. 14) of 0. 

All the other integral formula for Ji ~ /14 which are necessary in our calculation but not 
presented explicitly in the above can be taken from Appendix of . 

B. Anomalous dimension from each diagram 

In this appendix, we present the contribution to 7^^-* for hi from each Feynman diagram 
shown in Fig. 2, for which we use the notation 7^*^^ 7^^^ etc with obvious definition. We 
also use the symbols and etc to designate, respectively, the coefficients of the 

l/(167r2)2£:2 (double pole) term and the l/(167r2)25 (simple pole) term of the integrals Ji 
through J16 defined in the appendix of |^ and the appendix A above. 

As noted in the text, 7("'*'^'P) = 0, and ^(^'S.'^'''™-''?.'') are the same as those for fi calculated 
in 0]. Other results are as follows: {Cp = {N^ - l)/2N^, Co = A^c, Tr = Nf/2 for the color 
group SU{N,) ): 



7. 



16^ ^ 



(B.l) 



7^^) = "-^C,Ca. (B.2) 

Following the calculational procedure described in section 3, we obtained very cumber- 
some results for terms of the integral formulas, Ji ~ lie- It turned out, however, 
we can further simplify those results as was done in for fi. In the following, we present the 
simplified forms together with the expression with Ji ~ Jig. Although we had to calculate 
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the integral formulas for 7^£o separately from those for 7^>i, the final results for (i), (j), 
(k) can be written in a very simple form for all n, which supports the correctness of the 
calculation. 



^(0 

In 



16 (C| - ^CfCg) P-"- [p-^[ - I^in, 2) + 2/f^-(l, n + 1) + if ^'(n, 2) 
-27f ■^'•(1, n + 1) - 21^'' {n + 1,1) + 2/f ■^'•(n + 2, 0) - 21^'' {n + 2, 0) 
+27f -^(n, 2)} + p-'{inn, 1) - /f^(0, n + 1)- I^in, 1) + 27f ^(1, n) 

-1^(0, n+1)- 2/3^-^-(n, 1) + 2/3^-^> +1,0)} + 2/f^-(0, n) - 2/f/-(n, 0) 

+p-^{/f •''■(n, 1) - -^(O, n+1) - 27f -^(1, n) - I^in, 1) - -^(O, n + 1) 

+2/f -^ (l, n) + 4ir-{n + 1, 0)1 - /f -^'(n, 0) + /f '^-(O, n) + I^-''-{n, 0) 



2 n + 1 



Si{n) + 



n + 1 



(B.3) 



2p-2|/f ^-(n, 2) - 27f^(l, n + 1) - /f ^(n, 2) + 27f ^(1, n + 1) 
+2/3^-^-(n + 1, 1) - 2/f^-(n + 2, 0) + 2/|-^-(n + 2, 0) - 2/8^-^-(n, 2)| 
+p-^| - 37f -P(n, 1) + 37f ^(0, n + 1) + 3/2^-^(n, 1) - 47f ■^^(1, n) + /f ■^^ (0, n + 1) 
+4/3^-^-(n, 1) - 4/3^-^-(n + 1, 0)| + /f^-(n, 0) - 3/f^-(0, n) + il^f-{n, 0) 
+p-'[ - iF ^' in, 1) + 7f -^ (O, n + 1) + 4/f -^(l, n) + /f -^"(n, 1) + 37f -^(O, n + 1) 
-4/f -^(l,n) - 8/6^-^-(n+ 1,0)} + /f -^-(n,0) - 3/f-^(0,n) - /f ^■(n,0)" 



= 0. 

The contribution from (j) turned out to be identically zero for all n. 



(B.4) 



In 



8 Cf- oCfCg P 



p-\l + (-l)'^){/f-^ (n + 2, 0) - 27|-^ (n, 2) - 2/f/ (n, 2) 



+41- 



16 



n 



-I 



h 1) - 2/f/-(^ + 2, 0)| + p-'[inn. 1) + /f/-(n, 1) - /f/-(n + 1, 0) 



^f-{n, 1) + /r-(n + 1, 0) - /fr-(n, 1) + 4(1 + (-l)")/^!^, 1) 



-4(1 + (-l)")4^-(n + 1, 0)} + (1 + (-l)")/f-^-(n, 0) - (-l)"/f/-(n, 0) + 27f3^-(n, 0) 
+/f/-(n, 0) - 2(1 + (-l)")/f/-(n, 0) +p-'[inn, 1) + I^,^{n, 1) - I^/-{n + 1, 0) 



14 



-/fa ^-(n, 1) + iFs^'in +1,0)- I.Yin, 1)} - /f 0) + ^ (n, 0) + 0) 
_8(c^-^C.C.)(l + (-ir)(^-^ 



(B.5) 



This result is zero for odd-n. This can be easily proved by assigning the momenta as shown 
in Fig. 10. Then the denominator is symmetric under ki ^ k^. The numerator has a factor 
(^1 — ^2)"- It is convenient to consider the sum of Fig. 10 and its charge conjugation. Then 
one can easily show the contribution from (k) becomes for odd-n with the use of lemma in 
section 3. 



In 



n r f 

SCfCgP''' E P']- ^I'^'il -l,n-l + 2) + /f-^-(/, n-l + 1) 
1=1 '- 

+7f-^-(Z - 1, n - Z + 2) - 7f-^-(Z, n - Z + 1)| + 7f^-(Z, n - I) 

Sfin) S2{n) {2n-l)Si{n) 



= ACfCg 



S,{n) - 2S{n) 



2{n+l) 2{n+l) 



{n + iy 



+ 



2n 



(B.6) 



In 



= -16 (Cl - ^CfCg^ p-^ £ [ ^~'{ ~ + + I^in -1 + 1,1) 

+I^2''in -1,1 + 1)- 2I^f-{n -/ + !,/) + /f/-(n - / + 2, / - 1) 



■12 
tS.p., 



-I^s'^in -1,1 + 1) + 2^f\n -1 + 1,1)- I^f\n -1 + 2,1-1) 



+/f/-(n - /, / + 1) - /f/-(n - / + L /) 



lSf-^ri-l,l)-Iff\n-l + l,l-l) 



rS.P., 



-Iff in -1,1)- Itf {n -1 + 1,1-1) 



S.P., 



— —8 ( Cp — -CpCc 



4hSi(n) + 



+2Ss{n) + S'^ 



n 



2S,{n) ^A{l-i-m 



.2) n + 1 ' {n + lf 
In the above expression, we introduced the function 

S{n) = ±\^S,{3). 
j=i J 



- S2{n) 
2n 

{n + iy^ 



4:S{n) - 8S{n) 



(B.7) 



(B.8) 
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Figure captions 



Fig. 1 Basic vertices for • A with 0, 1 and 2 gluon lines. 

Fig. 2 Two-loop corrections to the one-particle irreducible two-point Green function which 
imbeds the operator ■ A. For diagrams not symmetric, crossed diagrams must be 
included. In Figs, (f) and (h), the contributions from the Faddev-Popov ghost and 
tadpole diagram must be added. 

Fig. 3 The same as Fig. 2 but for the flavor-singlet distribution. These diagrams become 
identically zero for hi. 

Fig. 4 Two-loop corrections to the quark self-energy. The ghost-loop and the tadpole contri- 
bution is implied as in Fig. 2. 

Fig. 5 Momentum labeling for the diagram Fig. 2(i). 

Fig. 6 The counter term contributions to Fig. 2(i). 

Fig. 7 The coefficients of the two-loop anomalous dimension ^l'^'^^^ for Nf — 3 (circle) and 
Nf — 5 (triangle). 

Fig. 8 (a) 7i{''^(i)/2/3i (circle) and 7^''*(°)/2/3o (square), (b) 7^''*(i)/2/3i - 7^''*W/2/3o for Nf = 
3, 5. 

Fig. 9 (a) The LO and NLO Q^-evolution of the tensor charge normahzed ai Q'^ — 1 GeV^. 
(b) The LO and NLO Q^-evolution of the first moment of hi{x,Q^) and fi{x,Q'^) 
normalized at = 1 GeV^. 

Fig. 10 Momentum labeling for the diagram Fig. 2(k). 
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